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Evidence for jet launching close to the black hole in 
GRB 101219B - a Fermi GRB dominated by thermal emission 

J. Larsson 1 , J. L. Racusin 2 and J. M. Burgess 1 


ABSTRACT 

We present observations by the Fermi Gamma-Ray Space Telescope Gamma-Ray Burst Mon¬ 
itor (GBM) of the nearby (z = 0.55) GRB 101219B. This burst is a long GRB, with an associated 
supernova and with a blackbody (BB) component detected in the early afterglow observed by 
the Swift X-ray Telescope (XRT). Here we show that the prompt gamma-ray emission has a 
BB spectrum, making this the second such burst observed by Fermi GBM. The properties of 
the BB, together with the redshift and our estimate of the radiative efficiency, makes it possible 
to calculate the absolute values of the properties of the outflow. We obtain an initial Lorentz 
factor T = 138 ± 8, a photospheric radius r p hot = 4.4 ± 1.9 x 10 11 cm and a launch radius 
ro = 2.7 ± 1.6 x 10 7 cm. The latter value is close to the black hole and suggests that the jet 
has a relatively unobstructed path through the star. There is no smooth connection between the 
BB components seen by GBM and XRT, ruling out the scenario that the late emission is due 
to high-latitude effects. In the interpretation that the XRT BB is prompt emission due to late 
central engine activity, the jet either has to be very wide or have a clumpy structure where the 
emission originates from a small patch. Other explanations for this component, such as emission 
from a cocoon surrounding the jet, are also possible. 

Subject headings: gamma-ray burst: general - gamma-ray burst: individual (GRB 101219B) - radiation 
mechanisms: thermal 


1. INTRODUCTION 

Insight into the physical properties of jets in 
gamma-ray bursts (GRBs) is dependent on identi¬ 
fying the emission mechanism responsible for the 
prompt gamma-ray emission. This is an unsolved 
problem, where no theoretical model has yet man¬ 
aged to explain the full range of observed param¬ 
eters for the smoothly broken power-law shape of 
the spectra. In particular, it is well known that 
synchrotron emission cannot e xplain the hard low- 


energy slopes of m any bursts (|Crider et al.lll997 : 
IPreece et al.l Il998h . while the majority of GRB 
spectra are significantly wider than expected for 
the simplest cases of thermal emission from the 
jet photosphere. From observations with Fermi 
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there is growing evidence that emission from the 
jet photosphere is important in at least some 
bursts, although it is usually seen to c o-exist with 


a dominant non-the r mal c o mponent (IRvde et al 


2010: Guiriec et al. 2011 


Guiriec et al. 12013 : iBurgess et al. 2014b). 


Axelsson et al.1 20121 


Due to the difficulties in discriminating between 
different emission models, observations of GRB 
spectra which are close to pure Planck functions 
are very important. There have previously been a 
few G RBs with such spectra observed by BATSE 
(jRvdel 2004 ) as well as one observed by Fermi 
GBM ( Ghirlanda et al. 20131 ). For these bursts, 
however, the redshifts and radiative efficiencies 
are unknown, leading to uncertainties in the ab¬ 
solute values of outflow parameters derived from 
the spectral fits. Here we present observations of 
GRB 101219B, the first GRB with a spectrum that 
is well described by a Planck function that also has 
a detected afterglow and a known redshift. 
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The prompt emission of GRB 101219B was 
observed by both Swift Burst Alert Tel escope 
(BAT) and Fermi G BM ( Gummings et al. 20101: 
van der Horstl l2010h , and the afterglow was de¬ 
tected from X -rays to optical energies by several 
observatories (Gelbord et al.1 2010t Olivares et al 


20ld : de TJgarte Postigo et al.l 201 lib The" opti¬ 
cal signature of an acco mpanyi ng S upernova (SN) 
was also detected dde TJgarte Postigo et al. 2011 


Olivares et al. 201 ill and the redshift was deter- 

mi ned to be 0.55 . An analysis of the SN spectrum 
by ISoarre et al. (feoill ) showed it to be a broad¬ 
line type Ic, with properties simil ar to SN 1998bw, 
the archetypal GRB-SN ( Galama et al . ITill . In 
addition, the Swift XRT spectra of the early af¬ 
terglow revealed the presence of a highl y signifi¬ 
cant BB component ( Starlin g et ah 20121) . Such a 
component has also been found in a small num¬ 
ber of other GRBs and the origin has been at¬ 
tributed variously to the SN shock breakout, emis¬ 
sion from the cocoon surrounding the jet, or late 
emission from the jet itself (IStarling et ah _ 2012 ; 
Snarre fc Starlingll2012t iFriis fc Watson 201 .ill . 

This paper is organized as follows: we present 
the spectral analysis of the gamma-ray observa¬ 
tions in section I2J1 estimate the radiative effi¬ 
ciency from the Swift XRT observations in section 
12.21 and derive the outflow parameters in section 
12.31 We discuss our results in section [3] Through¬ 
out this paper we assume a flat Universe with 
Hq = 71 km s" 1 and Om = 0.27. Error bars cor¬ 
respond to 1 sigma unless otherwise stated. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Gamma-ray spectral analysis 

In the analysis of the prompt phase of GRB 101219B 


we focus on the data from GBM (Meeganetal 


2009), which has a stronger signal as well as 
a wider energy range than the BAT. We do 
check the BAT data for consistency though. 
GRB 101219B is a relatively weak burst with 
a fluence of 3.99 ± 0.05 x 10~ 6 erg cm -2 in 
the 10 — 1000 keV energy ba nd and duration 
Xgo = 51.0 s ( von Kienlin et al. 2014 1. The light 
curve of the Nal detectors (8 — 1000 keV) is dom¬ 
inated by a single peak, as shown in Fig. [T] 

For the spectral analysis we use the Nal 3, 6 
and 7 detectors, which all have angles with re¬ 
spect to the burst of < 60° . There is no signal in 


the BGO detectors (200 keV - 40 MeV) or in the 
Large Area Telescope (LAT, 30 MeV - 300 GeV). 
The fits were performed using xspec vi2.8.i and 
pgstat statistics. To select time intervals for the 
spectral analysis we binned the light curve using 
the Bayesian Blo cks tech nique (as implemented in 


VanderPlas et al. 20141 ), which resulted in 


astroML, 

two time intervals, shown in Fig. IT] 

GRBs are typically well described b y the so- 
called Band function ( Band et al. 19931) . which is 
a smoothly broken power law characterized by a 
peak energy (A pea k) and power-law indices a and 
/3 below and above the peak, respectively. The 
results from fitting GRB 101219B with a Band 
function are listed in Table |T| The fits show that 
a is very hard, 0.63±0.05 and 0.88T(k26 i n the ^ rst 
and second interval, respectively, and close to the 
Rayleigh-Jeans value of 1. For comparison, the 
median value o f a f or long GRBs in the GBM cat¬ 
alog is —0.89 ( Gruber et al.ll2014 ) and synchrotron 
emission in the slow -coolin g regime has a = —2/ 3 
( Preece et al. 19981 see also lBurgess et al1l2014al) . 
In addition, this burst also has an unusually low 
Apeak (~ 70 keV in the first time interval), which is 
at the lowest 10th percentile of constra ined A pea k s 
for GBM bursts ( Gruber et al.l f2014h . and even 
more exceptional when considering the low red- 
shift. Motivated by these unusual spectral prop¬ 
erties we attempt to fit the data with a BB model, 
which we find provides a good fit to the spectra. 
The fit to the first time interval is shown in Fig. [2] 
and the best-fit parameters are presented in Table 
|T| The temperature of the BB is seen to decrease 
between the two time intervals. 

In order to test the hypothesis that the spec¬ 
trum of this burst, which has rather poor signal- 
to-noise, could still be consistent with slow-cooling 
synchrotron emission we use xspec to simulate 
10,000 spectra with a = —2/3 and /3 and A pea k 
taken from the fit to the first interval. The simu¬ 
lations are based on the background spectrum, re¬ 
sponse files and exposure time of the real data, and 
uses Poisson statistics to create synthetic spectra. 
Fitting these spectra with a Band function gives 
best-fit values of a > 0.63 in only six cases, corre¬ 
sponding to a probability of 6 x 10~ 4 of obtaining 
such a hard value of a by chance. The hypothesis 
that a = —2/3 can thus be rejected with very high 
confidence. 

Even for the case when the prompt emission 
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Fig. 1.— Light curve of GRB 101219B from Nal 3 
with 1 s bins. The two time intervals used for the 
spectral analysis are shown as dashed lines. The 
black line shows the background level, which was 
determined by fitting the light curve before —10 s 
and after 70 s with a straight line. 



Energy (keV) 

Fig. 2.— Fermi GBM spectrum of GRB 101219B 
from the first interval in Fig. |T] fitted with a BB 
model. Data from Nai 3, 6 and 7 are shown in 
black, red an blue, respectively. 


is completely dominated by the jet photosphere 
the observed spectrum is not expected to be a 
pure Planck function. Geometric effects are ex¬ 
pected to broaden the spectrum, with the hard¬ 
est a values of +0.4 —1-0.5 being obtained in the 
case of a wide jet viewed along the line of sight 
( Pe’ei 200ft : Beloborodov l201ll : iLundman et al 


20131 : Deng fc Zhanel 2014 1. In addition, we 


ex¬ 
pect some broadening due to time-evolution of 
the spectrum during the ~ 20 — 30 s time inter¬ 
vals from which the spectra were extracted (it is 
not possible to obtain good constraints on spec¬ 
tral parameters on a shorter time scale for this 
burst). A comparison of the Band and BB fits 
can be used to constrain the amount of broaden¬ 
ing in GRB 101219B. While a is close to the value 
expected for a pure BB, /3 = —2.60 is harder than 
the exponential cutoff of a BB. In line with this, 
a fit with a cutoff power-law model over-predicts 
the data at low energies as it needs to compen¬ 
sate for the somewhat harder high-energy slope. 
The flux derived from the BB fit to the first time 
interval is approximately 60% of the flux derived 
from the Band fit, which places an upper limit on 
the amount of broadening. As already mentioned, 
some of this will be accounted for by geometric 
effects and temporal evolution, but a small contri¬ 
bution from non-thermal emission, such as Comp- 
tonization of the BB, is also possible. 


For consistency we also analyzed the Swift BAT 
data of the burst. Spectra were created using 
standard BAT software (HEASOFT 6.15 and bat- 
GRBPRODuct). The results of the joint spectral fits 
are fully consistent with those presented above. 
We note, however, that the BAT spectra on their 
own are adequately fit with a simple power law 
(r = —1.7 + 0.2 in the first time interval). This is 
due to the combination of a weak signal and the 
limited band pass (15 — 150 keV), which means 
that the power law simply provides an average of 
the spectrum covering the spectral peak. 


2.2. Radiative efficiency and jet opening 
angle 

The radiative efficiency is defined as 77 = 
-E 7 i i so /(.E 7l iso + E\f), where E lt i so is the isotropic 
energy emitted in gamma rays and is the ki¬ 
netic energy of the jet. A measurement of 77 is 
needed in order to calculate the absolute val¬ 
ues of fireball parameters from the BB prop- 


3 
















































































































































































erties (see section m - The value of E k can 
be derived from X-ray after g low o bservations, 
as described by IZhang et al.l ( 2007 b We ob¬ 
tained the X-ray afterglow light curve and spectra 
for GRB 101219B from the Swift XRT reposi¬ 
tory ( Evans et al • §Q(S). Fitting the count rate 
l ight curve using the method of iRacusin et al. 
( 2009 ?) yields three power-law segments with 
temporal decay indices of a\ = 1.92 ± 0.07, 


«2 = 0.46lg 22 , «3 = 0.74^o'og with two breaks at 
tbi l = 4.1+i'g ks and R>r l = 581'gg ks. These seg¬ 


ments can be identified as late prompt emission, 
a shallow phase with continuous energy injection 
and the normal decay due to the decelerating fire¬ 
ball. No jet break is detected. A flare around 
300-600 s was modeled by a Gaussian. The light 
curve and the fit are shown in Fig. [3] 

We estimate E k from the normal decay phase, 
assuming typical values for the fraction of energy 
in electrons and magnetic fields (e e = 0.1 and et, = 
0.01) and neglecting Inverse Compton emission. 
We first determined the intrinsic column density, 
An.int, by fitting the full spectrum during the nor¬ 
mal decay (i.e. after 58 ks) with a power law modi¬ 
fied by intrinsic and Galactic absorption (with the 
Galactic absorption fixed at 3.1 xlO 20 cm -2 ). This 
yields Adpint = 1 . 2 + )' 2 x 10 21 cm -2 , consistent 
with the results of Starli ng et al. (2012). We then 
fit the model with Nh ^ n t fixed at this value to 
the spectrum at 58 ks and use the fit results to - 
gether with equation (5) in IRacusin et al.1 ( 2011 ) 
to derive E^. Assuming slow-cooling synchrotron 
emission with an electron spectral index p > 2 , 
the observed spectrum is consistent with being 
above the synchrotron cooling frequency. We ob¬ 
tain a value of E^ = 6.4 ± 3.5 x 10 52 ergs. With 
^ 7 , iso = 3.4 ± 0.2 x 10 51 ergs (over 10 — 1000 keV 
in the rest frame of the burst) this corresponds to 
a relatively low efficiency of r) = 5 ± 2%. 


This value of rj is associated with systematic 
uncertainties from our assumptions, in particular 
regarding the values of e e and £&. Although the 
values of e& reported in the literature span many 
orders of magnitude, we note that our estimate 
of Eb is only weakly dependent on this parame¬ 
ter. As an example, assuming a much lower value 
of eh = 5.5 x 10 ~ 6 (the median value found by 
Santana et al.1120141 for p = 2.8 and a constant 
density medium) gives 77 ~ 2 %, i.e. only slightly 
lower than our estimate above. The dependence 


on e e is stronger, but measurements of this pa¬ 
rameter have a rather narrow distribution, with 
~ 60% of GRBs fa lling in the range 0.1 — 0.3 
(ISantana et al. 2014^, in li ne w ith theoretical sim¬ 
ulations ( Sironi fe Spitkovskvl l201ll ). Increasing 
our value of e e to 0.3 at the high end of this dis¬ 
tribution gives 77 ~ 20 %, which gives an estimate 
of the uncertainty from this parameter. 


A limit on a jet break can be obtained from 
the last time when such a br eak could be presen t 
without being detectable (see lRacusin et all2009l) . 
For GRB 101219B this time is tb > 15.3 days. 
Assuming an ambient number density of 1 cm -3 
and 77 and F 7 j i so from above, this corresponds to a 
limit on the jet opening angle of 9j et > 17.1° and a 
collimation corrected energy E 7 > 1.5 x 10 5 ° ergs. 


2.3. Fireball parameters 


The observed temperature and flux of the BB, 
together with 77 and z, can be used to calculate the 
properties of the outflo w in the stand ard fireball 
scenario. Following IPe’er et al. ( 2007 ) we calcu¬ 
late the Lorentz factor (T), the radius where the 
fireball starts accelerating (ro), the saturation ra¬ 
dius (r s , where the acceleration stops) and the ra¬ 
dius of the photosphere (r p h). In the calculations 
we assume that there is no dissipation below the 
photosphere and that T is much larger than the in¬ 
verse of the opening angle of the jet. The results 
for the two intervals are presented in Table |T| 


The main source of uncertainty in these results 
is the radiative efficiency, as already discussed 
above. Another uncertainty is the contribution 
of non-thermal emission to the gamma-ray spec¬ 
trum. A very conservative limit on this can be 
obtained by comparing the fluxes from the Band 
and BB fits, as described in section 12.11 Taking 
the difference of these fluxes to be non-thermal 
emission would reduce the values of ro and r s by 
about a factor of two, while at the same time in¬ 
creasing the values of r p h and T by about 15%, as 
compared to the values presented in Table [l] 


3. DISCUSSION 

As noted in Section [T] there is growing evidence 
for photospheric emission in the prompt phase of 
GRBs, but very few examples where the spectra 
are close to pure BBs. It is clear that special condi¬ 
tions are required in order to obtain such spectra. 
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Table 1 

Spectral fits and Fireball parameters 


Time a 

(s) 

Band a b 

Band 0 h 

Band Ep eak 
(keV) 

BB kT c 
(keV) 

r d 

r 0 

(cm) 

r a 

(cm) 

r d 

r ph 

( cm) 

r d 

-5-18 

18-50 

0.63 ± 0.05 

0.88lS: 2 26 

—2 60 +O18 
Z.OU — o 37 

o qo + 0.12 

*- zo -0.22 

72.9 ± 14.8 
38.6+ 6 ' 2 

19 l+°- 7 

11 2"*"°' 8 
- L1 - z -0.7 

2.7 ± 1.6 X 10 7 
4.5 ± 2.8 X 10 7 

3.7 ± 2.3 X 10 9 
4.1 ± 2.6 X 10 9 

4.4 ± 1.9 x 10 11 
4.9 ± 2.3 X 10 11 

138 ± 8 

92 ± 8 


Note.— 

a Time interval relative to the GBM trigger. 

b Best-fit parameters from fitting the GBM data with a Band model. 
c Best-fit parameters from fitting the GBM data with a BB model. 
d Fireball parameters derived from the BB fits. 
u Parameter error range unconstrained. 



Fig. 3.— Swift XRT light curve of GRB 101219B. 
Cyan data points were taken in WT settling mode, 
blue points in WT mode and red points in PC 
mode. The solid black lines shows the fit with 
three power-law segments, separated by breaks 
marked by green dashed lines. The dotted black 
line shows a flare that was removed from the fit. 


First, a fairly wide jet viewed close to the line of 
sight is needed in order to avoid significant geo¬ 
metric broadening. This is in line with the limit 
on the jet opening angle from section 12.21 How¬ 
ever, it is interesting to note that the low-energy 
slope of GRB 101219B is harder than predicted for 
a spherically symmetric wind (see section mu), a 
possible explanation for this is that the emission 
arises from a small patch within a clumpy jet, in¬ 
validating the assumption of a steady flow that is 
spherically symmetric inside 1/T. An additional 
requirement for having a narrow spectrum is the 
lack of significant dissipation of energy in the jet. 

The outflow parameters in Table [l] provide im¬ 
portant information about the jet. It is partic¬ 
ularly interesting to relate the base of the jet to 
the size of the presumed compact remnant. In the 
first time interval rg = 2.7 ± 1.6 x 10 7 cm, which 
corresponds to ten times the event horizon for a 
18 ± 11 Mq maximally spinning black hole (or a 
9 ± 5 Mq non-spinning one). This mass range is 
close to the expected values for the generally as¬ 
sumed massive progenitor star, and suggests that 
the jet is launched from close to the black hole and 
then has an unobstructed path through the star. 
The latter requires a highly asymmetric explosion, 
in line with the observationally and theoretically 
well-established view that SN explosions are asym¬ 
metric. The only other GRB for which absolute 
values of the outflow parameters have been calcu¬ 
late^ is the highly energetic GRB 090902B, where 
Pe’er et all ( 2012( 1 found ?'o = 3.0 — 7.4 x 10 8 cm, 
i.e. significantly larger than for GRB 101219B. 
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That the conditions in this burst are different from 
GRB 101219B is also clear from the fact that the 
spectral shape is broader, most likely indicating 
subphotospheric dissipation. 


3.1. Connection with the blackbody com¬ 
ponent in the early afterglow 


GRB 101219B belongs to a small group of 
bursts in which a highly significant BB compo¬ 
nent has been det ected in Swift XRT data of 
the ea rly afteredow (IStarling et al.ll2012n . Specif¬ 
ically, IStarling et al. (2012 h fit'four time intervals 
between 180 and 1080 s after the GBM trigger 
(around the time of the flare in the light curve, 
see Fig. 0 using a model comprising a power law 
and a BB, modified by absorption (see their table 
5 for the full fit results). The BB has an initial 
rest-frame temperature of 0.3 keV and then de¬ 
creases with time. 


In order to investigate any possible connection 
between the BB in the prompt emission and the 
BB in the early afterglow we plot in Fig. [4] the 
evolution of the BB temperature as a function of 
time for the two datasets. The temporal evolution 
of the temperature in the Swift XRT data is well 
described by a power law (shown by the dashed 
line, with decay index —0.85), but the extrapola¬ 
tion of this fit to early times is clearly inconsistent 
with the GBM data. A broken power-law model 
also fails to connect the two data sets. In Fig. [I] 
we also plot the parameter 1Z = (-Fbb/ctT 4 ) 1 / 2 , 
which is related to the effective transverse size of 
the photosphere. This parameter is also seen to 
evolve differently during the two phases. 


The lack of a smooth connection between the 
observed properties of the BB in the two data sets, 
together with the fact that the BB flux a ctually 
increases at late times ( Starling et al. 1201 21. rules 
out the scenario that the late emission is due to 
high latitude emission from the jet after the cen¬ 
tral engine has died. We therefore explore the sce¬ 
nario that the late emission is due to late central 
engine activity, and calculate the outflow param¬ 
eters from the low-energy BB as described in sec¬ 
tion 12.31 The resulting values of T and r rp hot are 
shown in Fig. 0 for two scenarios: (i) all the late 
emission is prompt emission and (ii) the BB is the 
prompt emission from the jet, while the power-law 
component arises in a different region where the 
jet has started interacting with the circumstellar 


medium. The two scenarios produce comparable 
results. We see that at late times the jet has a low 
r in the range 3 — 6. In order for such a slow jet to 
still produce a narrow BB spectrum the jet has to 
be very wide (^ 20°, consistent with the limit of 
9j > 17.1° from section I2.2[) or, alternatively, the 
emission could originate from a small patch within 
a clumpy jet. 

As already discussed by Starlin g et ab (2012), 
there are also other possibilities for the origin 
of the low-energy BB. In particular, they find 
that emission from a cocoon surrounding the jet 
can explain the main features of the emission, 
while shock breakout from a SN is disfavored 
from the large inferred radius (~ 10 12 — 10 13 cm, 
much larger than the 10 11 cm expected from the 
presumed Wolf-Rayet progenitor). While break¬ 
out through a thick wind could accommodate 
the large radius in the SN scenario, the spec- 
trum is in this case e xpected to be flat in vF v 
( Svirski fe Nakar 2014f) . contrary to what is ob¬ 
served for GRB 101219B. 


This work was supported the Swedish National 

Space Board. 
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